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Radial wave functions for power-law potentials are approximated with the help of power-law sub- 
stitution and explicit summation of the leading constituent WKB series. Our approach reproduces 
the correct behavior of the wave functions at the origin, at the turning points and far away from 
the turning points 
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INTRODUCTION 



where 



Many interesting physical problems require solving the 
Schrodinger equation for spherically symmetric poten- 
tials V(r). The separated radial Schrodinger equation 
can be written in the form 



h 2 d 2 iP(r) 
2m dr 2 



+ (V(r) + V c (r)-E)ip(r)=0 (1) 



which is identical to the one-dimensional Schrodinger 
equation with an effective potential given by the sum 
of the origin potential V(r) and the centrifugal potential 



V c (r) 



h 2 l(l + 1) 

2mr 2 



(2) 



In the present work we consider power-law potentials 

V(r)=ar k , a > 0, k > 1 (3) 

which are very important in particle physics. It is known 
that only several potentials permit the exact solutions of 
the Schrodinger equation. Usually realistic physical cal- 
culations are impossible without different approximation 
techniques. One of the earliest and simplest methods of 
obtaining approximate eigenvalues and eigenfunctions of 
the radial Schrodinger equation is the WKB method (see, 
e -g-j [H>0 an d references therein). 

It is known [J, Q that a suitable transformation of 
the initial equation improves results of an approximation 
technique. We examine the power-law substitutions 



= q s , s> 0, ip(r) 



-^/ 2s ^(q). (4) 



The transformed equation is 



-h 



d 2 q{g) 
dq 2 



Q(q)*(q)=0 



(5) 
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Q(q) 



2ms 

h 2 



aq 



(2+k)s 



Eq 



2 s 



\s\l + l/2Y-~ . 



(6) 



Of course, the exact solutions do not depend on some 
substitution. However, we are interesting in the approx- 
imate solutions. R. E. Langer [l[ has used the particular 
case s = 2 when he was applying the WKB approxima- 
tion to the Kepler problem. 

The WKB approach deals with the logarithmic deriva- 
tive 



Y(q) 



dln*(g) 
dq 



(7) 



that satisfies the nonlinear Riccati equation 

-h 2 {^ + Y 2 {q)^ + Q(q)=0 (8) 

where Q(q) is an arbitrary function of q (naturally not 
only the special expression (6)). The WKB series 



±Ql/2 + y>F± (g ) 



(9) 



are the asymptotic expansions in powers of Plank's con- 
stant H of two independent particular solutions of the 
Riccati equation. The usual WKB approximations 



Y^ KB (q) = h- 1 \ ±QV*+J2h n Y±( q )\ 

contain a finite number of leading terms YJp (q) from the 
complete expansions Y^(q). These approximations are 
not valid at the turning points where Q(q) = and at the 
origin q — 0. While in most cases of improvements of the 
WKB method (see, e.g., 0y,13;13l) the main purpose is 
to achieve highest accuracy in eigenvalue calculation for 
the radial Schrodinger equation, our goal is to construct 
satisfactory approximate eigenfunctions with the correct 
behavior at the origin, at the turning points and far away 
from the turning points. 
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II. NEW APPROXIMATIONS FOR 
LOGARITHMIC DERIVATIVES 

The analysis of the well-known structure of the leading 
Y^ (q) and recursion relations [1, 0] allows us to recon- 
struct the asymptotic WKB series as the infinite sums 



Y±(q) = ±h- 1 Q 1 / 2 



(10) 



of new constituent (partial) asymptotic series Z^ s 7 - (<?) 

The complete series (q) are approximated by a finite 
number of leading constituent series Z^ s j(q) in contrast 
to the use of a finite number of leading terms Y^ 1 (q) in 
the conventional WKB approach. Using notation 



a(q) 



Q(q) 



K 2 ' 3 IdQ^/dq] 2 / 3 ' 



(11) 



dVa S , 2 , „. ± ± _ 1 (r, d yt.l 



da ^^ = ~A 2a -da 



y«s 



(16) 



Eq.(15) is the Riccati equation for the logarithmic deriva- 
tives of linear combinations of the well-studied Airy func- 
tions Ai(a) and Bi(a) [l(| • 

The particular expressions 



yf(a) = ^ln (Ai(o)±iBi(o)), 



J/i (a) = ^lnAi(a), yf(a) = -^-hiBi(a) 



h(q) = 



dQ(q)/dq 



H 2 / 3 \dQ(q)/dq\ 2 / 3 ' 



(12) 



correspond to the conventional WKB series . 
Generalization consists in the use of formulas 



b 2 (q) 



d 2 Q(q)/dq 2 
dQ(q)/dq 



(13) 



yi (a; t) = — In (Ai(a) + tBi(a)) . 
da 



(17) 



we are able to rewrite two first leading constituent series 
in the form 



±h- 1 Q^ 2 + Zt sA (q)+Zt s ^q) 
= h(q)yt,i( a ) + b z(<l)ytsA a )- 



(14) 



Direct verification shows that the series yf s x (a) and 
Vas 2 ( a ) satisfy equations 



V2(a;t) = — [~8a 2 ( yi (a;t)) 2 - 3 - 4a yi {a;t) + 8a 3 ] 

(18) 

with a mixture parameter t. As a result we get new 
approximate logarithmic derivative 



Y avp {q)=Y(q-t) = b l (q)y 1 {a;t) + b 2 (q)y 2 {a;t) (19) 



dyt, 



da 



- + d£,i) a = 



(15) which satisfies the following equation 



h 2 (^M. • Y 2 (q;t)) QUP 



+ 



V dq 
d 2 Q(q)/dq< 
dQ(q)/dq 



( d 3 Q(q)/dq 3 
V dQ(q)/dq 



V2{a\t) 



8 4 1 

vl{a;t) - -y 2 (a;t) + -a yi (a;t)y 2 {a;t) - - 



(20) 



instead of Eq. (8). 



tion near the turning points. 



It is not surprising that the asymptotics of our approx- 
imation coincide with the WKB asymptotics far away 
from the turning points. At the same time our approxi- 
mation reproduces the known 0] satisfactory approxima- 
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III. NEW APPROXIMATIONS FOR WAVE 
FUNCTIONS 

Now we can construct the approximate radial wave 
functions for the bound states in the case of a power- 
law potential when Q(q) is of the form (6). 

First, we must reproduce the correct limiting behavior 
at the origin. In this case we have the following exact 
expressions (r — > 0, q — > 0) 



ipex(r) 



„si+(s+l)/2 



Y ex (q) 



sl + {s + l)/2 



(21) 



At the same time we can derive relations 

/s 2 1/3 
a(q)^a = f— (i + 1/2) 2 - — 



Y(q;t) -» --(2a m (a ;t) + 3y 2 (a ;t)) (22) 

in the framework of our approach. We obtain the alge- 
braic equation for determining the value of t. Its solution 
is 



to 



—c(l, s)Ai(ao) + a (dAi(a ) / 'da ) 
c(l, s)Bi(ao) — ao(c?Bi(ao)/(iao) 



(23) 



where 



c(l,s) = 1 




+ s(l + 1/2) + 1/2 



Two real turning points g_ and q+ (Q(q±) = 0) sepa- 
rate three regions. 

In the first region where < q < q- we select 
the unique approximate particular logarithmic deriva- 
tive Y(q;to). In the second region where q- < q < q + 
we must describe the oscillatory solution of the original 
Schrodinger equation (1). Therefore in this case we se- 
lect two approximate particular logarithmic derivatives 
Y(q; +i) and Y(q; —i). In the third region where q > q+ 
we must describe only the decreasing solution of the orig- 
inal Schrodinger equation (1). Therefore in this case 
we select the unique approximate particular logarithmic 
derivative Y(q; 0). Note that in the case ( = 0,s = lwc 
put q- = 0. 

Since the turning points arc ordinary nonsingular 
points in our approach, no question of connection formu- 
las arises in contrast with the conventional WKB method. 
Matching particular solutions at the turning points we 
obtain the continuous approximate radial wave function 



ipap P {r) = N app A s ~ 1)/2s ^ app (q) 



(24) 



where ^ ap p(q) is represented by the following formulas 
= cos0exp J* Y(q';to)dq'^j (25) 

if < q < q- , 



if g_ < q < q + , 

•.(.) - i(-D"-P (£ r«.,ow) exp ^ m±ii±mzii ^ 

if q > q+ ■ , 



Here q = r 1 ^,^ — § — arctani , e — 
(dQ(q)/dq) \dQ(q)/dq\^ 1 and N app is a normalization 
constant. 

We have the new quantization condition 




Y(q;+i)-Y(q;-i) 



= ir(n + -) + cf>, n = 0,1,2... (28) 



which determines the spectral value of E implicitly. 

Note that up to now a value of s is not fixed. Numerical 
experiment for the power-law potentials (k > 1) shows 
that the best choice for I = is s = 1 and the satisfactory 
common choice for all Z > is s = 2. Thus the approx- 
imate eigenfunctions are determined completely and we 
can perform verification. 
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IV. APPLICATION TO POWER-LAW 
POTENTIALS 

It is convenient to test our approximation with intro- 
ducing the dimensionless quantities 

/2ma\ 1/{k+2) 



2 m 

h 2 a 2 / k 



k(/k+2) 



E. 



Then the Schrodinger equation is rewritten in the form 

Hip(x) - e<ip(x) = 



with the Hamiltonian 



H = - 



dx 2 



(29) 



First, we estimate our approximation in the case of the 
harmonic oscillator potential V(r) — ar 2 for which the 
exact wave functions ip ex (x) are well known Fig- 
ures 1, 2 and 3 show that the proposed approximation 
gives fairly accurate wave functions in this case. Here 
solid lines reproduce Tp a pp(x) and dashed lines reproduce 
W(ipapp(x) — ipex(x))- Note that without factor 10 the 
difference (ip ap p(x) — ^ ex (x)) is invisible in comparison 

with Ipapp(x). 

We can calculate the expectation values 



and 



a 



(e'app) 2 =< ^app\H 2 \ll), 



> 



J app 



> 



(30) 



(31) 



with the help of the normalized approximate wave func- 
tions (< ip a pp\ipapp >= 1)- It should be stressed that 
{zapp)) 2 7^ ( e 'app) 2 when the wave functions are not ex- 
act. Now we define the relative discrepancy 



d = 



l. 



We also calculate relative virial error 



^ Ipapp | 



dx' 2 



app 



> 



> 



- 1 



(32) 



(33) 



J app | 2 I Yapp 

which is equal to zero for the exact solutions. Finally, 
we characterize our approximation by the usual relative 
energy error 



Se = 



1 



(34) 



where e ex is the exact energy value. 




FIG. 1: Harmonic oscillator radial wave function 
for I = 0, n = 0. 




FIG. 2: Harmonic oscillator radial wave function 
for I = 1, n = 1. 




FIG. 3: Harmonic oscillator radial wave function 
for I = 2, n = 2. 



We calculate the values Se for the linear and quartic po- 
tentials with respect to the accurate numerical eigenener- 
gies [lH- We compare our results with the recent results 
in the framework of an integral semiclassical method for 
calculating the spectra for spherically symmetric poten- 
tials [lH]. In the case of linear potential with / = our 
approach gives exact eigenfunctions and eigenenergies. 
Table 1 demonstrates validity of our approximation in 
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TABLE I: Numerical verification of the proposed approximation in the case of power-law potentials V(r) = ar k . 
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the linear, quadratic and quartic cases . 

V. CONCLUSION 

An old problem in semiclassical analysis is the devel- 
opment of global uniform approximations for the wave 



functions. In the present paper, this problem has been 
solved by means of the reconstruction of the WKB series 
and subsequent explicit summation of the leading con- 
stituent (partial) series. Such approach yields satisfac- 
tory description of the wave functions in the important 
case of the radial Schrodingcr equation with the power- 
law potentials. 
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